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ABSTRACT 


We  consider  here  a model  of  traffic  flow  on  a road  network. 

For  each  ordered  pair  of  nodes  there  is  a demand  function  which 
expresses  travel  demand  between  the  two  nodes  as  a function  of  travel 
times  on  the  network.  Each  road  (arc)  has  a delay  function  vjhich 
expresses  travel  time  on  that  arc  as  a function  of  total  traffic  flow. 
Our  objective  is  to  show  how  an  equilibrium  of  travel  times,  flows, 
and  demands  may  be  computed  under  conditions  which  are  simple,  general, 
and  plausible. 

To  solve  the  netvrork  problems  we  develop  techniques  for  solving 
the  stationary  point  problem.  These  techniques  for  the  stationary 
point  problem  are  the  best  of  which  we  are  avrare. 
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CHAPTER  1 


INTRODUCTION 

Let  ,&)  be  a directed  network  with  nodes  i in  9?  S'^d 
arcs  ij  in  tf.  Our  interest  is  in  the  traffic  flows  and  travel  times 
on  the  network. 

For  each  arc  ij  we  are  given  a delay  function  f . . which 
expresses  travel  time  on  arc  ij  as  a function  of  the  traffic  flows  on 
the  arcs  of  the  network.  We  might  expect  that  travel  time  along  arc 
ij  TTOuld  depend  on  the  flow  on  that  arc,  but  it  may  depend  on  flows 
along  other  arcs  as  well.  In  particular  a two-way  street  could  be 
modeled  as  a pair  of  opposing  arcs,  and  the  flow  on  one  of  the  arcs 
might  affect  the  travel  time  on  the  other  arc. 

Another  way  in  which  the  interactions  between  different  arcs 
may  arise  is  in  the  modelling  of  an  intersection.  Florian  and  Nguyen 
[10  1 employed  two  methods  of  representing  an  intersection,  see 
Figures  1 and  2.  Travel  time  across  the  intersection  v/ill  depend  on 
the  traffic  in  other  directions.  For  instance  the  time  required  to 
make  a left  turn  will  be  affected  by  the  density  of  oncoming  traffic. 

For  each  pair  of  distinct  nodes  i,k  we  are  given  a travel 

demand  function  g,  , which  expresses  demand  for  travel  from  i to  k 

1,  K 

as  a function  of  travel  times  between  nodes  on  the  network.  Demand  for 
travel  from  i to  k would,  of  course,  depend  on  travel  time  from 

I 


I 


Figure  1 


2 


i to  k.  But  it  might  also  depend  on  travel  times  between  other  pairs 
of  nodes,  for  instance  from  i to  some  alternate  destination.  Later  in 
this  chapter  we  will  address  some  of  the  issues  involved  in  deriving 
these  demand  functions. 

The  ability  to  hajidle  these  interactions  between  various  arcs 
and  origin-destination  pairs  is  the  principal  contribution  of  this  paper. 

Many  models  very  much  like  the  model  presented  here  have  been 
developed  over  the  last  25  years.  The  earliest  examples  of  this  kind 
are  due  to  Wardrop  [13]  and  Beckman  [2].  These  models  are  based  on  the 
notion  that  the  road  system  has  a large  community  of  users,  with  each 
user  tailing  the  quickest  route  available,  given  the  actions  of  other 
users.  In  addition,  the  number  of  trips  made  may  depend  on  the  time 
required  to  msike  a trip.  On  the  other  hand,  travel  time  on  a particular 
road  will  depend  on  traffic  volume. 

Numerous  solution  procedures  have  been  proposed  for  these  models. 
The  usual  approach  has  been  to  reformulate  the  equilibrium  problem  as  a 
convex  programming  problem.  This  has  the  advantage  that  convex  programming 
techniques  when  applied  to  this  type  of  problem  can  be  very  efficient  [9]. 
Rather  strong  conditions  must  be  satisfied  by  the  functions  f and  g in 
order  for  these  procedures  to  be  valid.  In  particular,  f must  be 
integrable  and  g must  have  an  inverse  which  is  integrable.  This  will 
be  true  if  each  f^^  depends  only  on  the  total  flow  y^^^  ^ and  each 
g depends  only  on  t and  is  monotone  decreasing.  The  equilibrium 
problem  is  solved  under  conditions  like  these  by  Beckman  [5],  Dafermos 
and  Sparrow  [5],  Florian  and  Nguyen  [9]>  and  others. 


5 


Dafennos  discusses  a model  in  which  travel  times  on  one  arc 
may  depend  on  the  flow  on  another  arc  in  order  to  model  the  interaction 
of  traffic  going  in  opposite  directions  on  a two-way  street.  In  order 
to  use  the  nonlinear  programming  approach  she  assumes  that  the  delay 
functions  are  continuously  differentiable,  with  symmetric  Jacobian 
matrices.  These  assumptions  are  essentially  a restatement  of  the 
integrability  conditions  mentioned  above.  These  are  reasonable  when  the 
street  is  s.ymmetric  and  the  traffic  flow  is  about  the  same  in  each 
direction.  Otherwise,  the  methods  presented  here  would  seem  more 
applicable.  This  would  be  particularly  true  for  the  modeling  of  an 
intersection  where  the  situation  would  normally  be  quite  unsymmetric. 

In  practice  the  demand  and  delay  functions  f and  g are  at 
best  empirical  fits  and  can  be  endowed  with  these  or  any  other  restric- 
tions which  may  seem  useful.  It  should  be  noted,  however,  that  the 
approach  used  here  does  not  require  such  restrictions,  and  these  restric- 
tions can  be  expected  to  complicate  the  handling  of  interactions  between 
various  arcs  and  origin-destination  pairs  as  described  above.  Indeed  to 
compute  an  equilibrium  it  is  only  necessary  that  the  delay  functions 

f. . be  positive  on  each  arc,  that  the  travel  demand  functions  g. 
ij  i>K 

be  nonnegative  and  bounded  for  each  pair  of  nodes,  and  that  the  network 
be  complete  that  is,  there  is  a directed  path  from  every  node  to  every 
other  node.  Unfortunately  we  must  pay  a price  in  computational  efficiency 
for  the  additional  generality.  For  problems  which  can  be  solved  by 
convex  programming  techniaues  that  would  be  the  best  approach. 

Next  we  present  the  mathematical  conditions  for  a user 
ecuilibrium.  The  travel  time  from  node  i to  node  k will  be  written 


'i,k- 


The  flow  on  arc  ij  with  destination  k will  be  written  y 


We  say  that  the  vectors  of  travel  times  t and  flows  x and  y are 
in  eauilibrium  if  the  following  conditions  hold. 


l.a) 

i k. 

i,  ke?? 

l.b) 

y > 0 

l.c) 

t.  <f..(x)  + t. 

' j,k 

V = 0 

k,  k 

i 5/  k. 

ij  £ a,  k € ^ 

l.d) 

^ij,k^^ij^*^  ^j,k  “ ^i,k^  " ® 

i 5^  k. 

ij  - a,  k e 7^ 

l.e) 

ij  € a 

It  may  be  useful  to  think  of  this 

system  as  a 

multicommodity 

network,  where  all  of  the  traffic  destined  for  a particular  node  k 
is  a separate  commodity,  all  of  which  must  be  shipped  to  node  k via 
the  network.  In  this  way  g,  ^(t)  is  the  amount  of  commodity  k 
which  must  travel  from  node  i to  node  k.  This  trip  will  traverse  a path 
of  arcs  from  i to  k. 

Condition  (l.a)  is  the  conservation  of  flow  equation.  It  says 
that  the  traffic  leaving  node  i with  destination  k is  the  sum  of  the 
traffic  arriving  at  node  i with  destination  k and  the  traffic 
originating  at  i with  destination  k. 

Condition  (l.b)  says  that  traffic  flows  cannot  be 
negative. 
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Conditions  (l.c)  and  (l.d)  require  that  traffic  flow  by  the 


fastest  route  available.  In  (l.c)  we  . jquire  that  t.  not  exceed 

i,k 

the  minimum  travel  time  from  i to  k given  the  flows  y on  the  net- 
work. Then  (l.d)  says  that  traffic  may  flow  only  along  a route  which 
achieves  this  minimum  travel  time.  Together  (l.c)  and  (l.d)  imply 
the  principle  of  minimum  travel  time.  This  says  that  if  any  traffic 
flows  from  i to  k,  that  is  if  Z.  y,  . ^ > 0,  then 


t 


i,k 


min  (f. .(y)  + t > 

A J > ^ 


Equation  (l.e)  relates  the  basic  flows  y to  the  total  arc  flows  x. 

More  generally,  the  functions  f and  g may  be  point-to-set 
maps.  In  this  case  we  say  that  travel  times  t,  flows  x and  y,  and 
arc  delays  t are  in  equilibrium  if  the  following  conditions  hold. 


J u 


l.b)  y > 0 

l.c)  T £ f(x) 


t.  , < T.  + t . 
i,k  - ij  3,k 


"kk 


yij,k(^ij  ^ ^j,k  - \k) 


l.e)  x^^  - Z 


i / k,  i,k  £ ^ 


i 9^'  k,  i j £ a,  k € ^ 

i 5^  k,  ij  £ 8,  k £ . 
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In  the  following  theorems  we  require  that  the  delays  and  demands 
be  expressed  in  terms  of  upper  semicontinuous  convex  valued  point-to-set 
maps.  In  practice  these  firnctions  are  usually  continuous  and  hence 
satisfy  the  requirements. 

THEOREM  1.  Suppose  is  a complete  network.  Suppose 

i)  f,  the  delay  function  is  positive,  convex  valued,  and  upper  semi- 
continuous  on  {x|x  > 0}  and 

ii)  g,  the  demand  function  is  nonnegative,  bounded,  convex  valued,  and 
upper  semicontinuous  on  {t|t  > 0}. 

Then  a solution  to  the  equilibri\uin  problem  exists  and  ran  be  computed  by 
the  Eaves-Saigal  algorithm.  □ 

Recently  Aashtiani  [1]  has  shown  how  the  equilibrium  problem  can 
be  formulated  as  a complementarity  problem.  Complementarity  problems 
will  be  defined  in  Section  1 of  Chapter  2.  Aashtiani' s approach  is  like 
the  approach  used  here.  He  proved  an  existence  theorem  which  is  similar 
to  Theorem  1,  but  he  shows  that  his  solution  procedure  may  not  converge 
under  the  conditions  of  his  theorem. 

The  conditions  of  Theorem  1 would  seem  to  be  sufficiently  general 
to  allow  the  modelling  of  a traffic  network.  However,  many  actual 
systems  have  been  modelled  with  demand  functions  which  are  unbounded  when 
the  corresponding  travel  times  are  small.  Beckmann  [3]  suggests  the 
following  formula  for  travel  demand. 
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®i  " ''■Pi'  : i: 


Here  p is  the  vector  of  populations  at  various  nodes,  and  a and  c 
are  given  constants.  Utility  functions  can  also  be  used  to  derive 
travel  demands.  In  neither  case  are  the  conditions  of  Theorem  1 satisfied. 
It  will  be  easy  to  see  that  the  explicit  function  given  above  satisfies 
the  conditions  of  Theorem  2.  Theorem  3 ar.d  the  discussion  preceding 
it  show  how  a demand  function  may  be  derived  using  utility  functions 
which  satisfies  the  conditions  of  Theorem  2. 

Here  t > € means  no  component  of  t is  less  than  e. 


THEOREM  2.  Suppose  is  a complete  network.  Let  e > 0,  Suppose 

i';  f,  the  delay  function  is  upper  semicontinuous  and  convex  valued, 

ii)  Each  f. .(x)  > e for  all  x > 0,  and 
ij  - 

iii)  g,  the  demand  function  is  nonnegative,  bounded,  convex  valued  and 
upper  semicontinuous,  on  {t|t  > e1 . 

Then  a solution  to  the  equilibrium  problem  exists  and  can  be  con^juted  by 
the  Eaves-Saigal  algorithm.  □ 


Theorem  1 is  just  the  special  case  of  Theorem  2 -vhere  e = 0. 
The  following  paragraphs  will  show  how  the  assumptions  on  g in 
Theorem  2 arise  naturally  when  g is  derived  as  the  sum  of  demands 
by  consumers,  each  maximizing  his  own  utility  function. 


Q 


We  now  suppose  that  the  network  is  used  by  several  consumers  j. 

Each  obseirves  the  condition  of  the  network  as  measured  by  the  travel 

time  vector  t.  Each  sets  up  a new  potential  demand  vector  that  maximizes 

his  utility.  His  utility  is  affected  by  the  amount  of  travel  between 

the  various  node  pairs  on  the  network,  and  the  total  time  devoted  to 

travel.  This  utility  is  expressed  through  a function  For  simplicity 

we  assume  that  each  consumer  has  an  upper  bound  on  the  total  amount 

of  time  he  has  available  for  travel.  The  sum  over  all  j of  the  demands 

for  travel  originating  at  i and  destined  for  k represents  g.  ^(t'. 

1,  k 

Given  travel  times  the  vector  of  t,  the  jth  consumer  will  find 
a pair  which  maximizes 

subject  to  ■jT'^-t  < 

> 0,  q^  < q^ 

where  is  the  vector  of  travel  demands  for  various  node  pairs,  and 

is  the  total  amount  of  time  devoted  to  travel  by  j.  The  vector 
product  v^-t  is  the  sum  of  the  time  spent  in  travel  between  all  node 
pairs,  hence  its  terms  represent  an  allocation  of  q^. 

The  following  theorem  shows  that  the  demand  function  derived  as 
the  sum.  of  these  individual  consumer  demands  satisfies  the  requirements 
of  Theorem  2.  It  tells  us  that  if  g is  derived  as  above  and  f > € 
for  any  f > 0 then  the  hypotheses  of  Theorem  2 will  be  satisfied. 
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THEOREM  Suppose  that  each  is  convex  and  continuous.  For  each 

t > 0 define  g(t)  = Z.  g'^(t).  Then  for  each  e > 0,  g is  non- 

J 

negative,  bounded,  convex  valued,  and  upper  semicontinaous  on 
{t|t  j t K □ 

In  Theorems  1 and  2 conditions  (l.a)  through  (l.e)  were  shown 
to  have  a solution  and  that  the  Eaves-Saigal  algorithm  can  construct 
a solution.  For  the  special  cases  mentioned  above  where  the  problem 
can  be  solved  as  a convex  program  the  equilibrium  solution  can  also  be 
shown  to  be  unique.  The  following  theorem  extends  the  uniqueness 
property  to  a more  general  situation. 

A point-to-set  mapping  f is  called  monotone  on  K C if 

12  11 
for  each  pair  of  vectors  x and  x in  K and  every  y C f(x 

^ 2 12  12 

ana  y'  £ f(x  ),  (x  -x  ) • (y  -y  ) >0.  The  mapping  f strictly  monotone 

1 2 

if  for  every  pair  of  distinct  vectors  x and  x in  K,  and  every 

y^  £ f(x^)  and  y^  £ f(x^),  (x^-x^) • (y^-y^)  > 0. 

The  following  theorem  includes  the  convex  programming  unique- 
ness result  as  a special  case  since  the  gradient  of  a strictly  convex 
function  is  strictly  monotone. 

THEOREM  h.  Suppose  in  addition  to  the  hypotheses  of  Theorem  2 that  f 
and  -g  are  strictly  monotone.  Then  the  set  of  solutions  to  the  equi- 
librium problem  is  convex  and  the  total  flows  x are  unique.  □ 

The  proofs  of  the  theorems  of  this  chapter  will  be  delayed  until 
Chapter  5 since  they  require  the  background  material  of  Chapter  2. 
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CHAPTER  2 


THE  STATIONARY  POINT  PROBLEM 

The  stationary  point  prohlem  provides  a common  framework  which 
includes  the  linear  and  nonlinear  complementarity  problems  as  well  as 
fixed  point  problems.  Applications  come  from  such  diverse  areas  as 
economics,  game  theory,  mathematical  programming  and  mechanics.  For  this 
reason  the  results  of  this  chapter  should  be  of  interest  beyond  the 
specific  equilibrium  problems  of  this  dissertation.  However  the  approach 
used  here  is  largely  dictated  by  the  problems  at  hand. 

In  Section  1 of  this  chapter  after  giving  the  necessary  basic 
definitions,  we  explore  some  of  the  more  important  special  cases  of  the 
stationary  point  problem.  Section  2 is  devoted  to  existence  theorems; 
in  Section  3 we  present  two  theorems  which  address  the  question  of 
uniqueness  of  stationary  points. 

2.1.  Definitions  and  Examples 

Let  D be  a subset  of  R*^.  A point-to-set  map  f from  D to 
subsets  of  R^  is  called  upper  semicontlnuous  if  x^  ->  x and  y^  2 f (x^) 
imply  that  y^  has  a cluster  point  in  f (x) . (We  must  make  a slightly 
stronger  definition  here  than  is  usual  since  we  do  not  assume  that  the 
image  set  is  compact. ) 

We  say  that  f is  lower  s eml cont i nuous  if  x^  -+  x and  y £ f (x) 
imply  that  there  is  a sequence  y^  -►  y such  that  y^  2 f(x^)  for  all  i. 
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If  f is  both  upper  semicontinuous  and  lower  semicontinuous  then  f 
is  continuous . 

The  graph  of  f is  {(x,y)|y  1 f(x)].  A necessary  and  sufficient 
condition  for  f to  be  upper  semicontinuous  is  that  the  graph  of  f is 
closed  and  f is  bounded  on  compact  subsets. 

If  f(x)  is  any  single  point  for  each  x then  the  above  definition 
of  upper  semicontinuity  coincides  with  the  usual  definition  of  continuity. 
The  same  is  true  for  lower  semicontinuity. 

If  D is  a subset  of  then  D*  is  the  set  of  all  nonempty 

compact  convex  subsets  of  D. 

n n* 

Let  K be  a convex  subset  of  R . Let  f:K  ->  R be  upper  serai- 
continuous.  A pair  (x,y)  is  a stationary  point  of  (f, K)  if  x £ K, 
y £ f(x),  and  (u-x)*y  > 0 for  all  u 6 K. 

As  an  example  consider  the  programming  problem 

min  g(x) 
xSK 


where  g is  differentiable. 

We  now  show  that  for  every  solution  x*  to  this  minimization 
problem  that  (x*,Vg(x*))  is  a stationary  point  of  the  pair  (7g,K), 
Suppose  X*  solves  the  minimization  problem.  Then  x*  € K.  Let 
u £ K.  Singe  g is  differentiable  we  may  expand  it  in  a Taylor  series 
for  A £ 1] . 

g(Xu  + (l-X)x*)  = g(x*  + X(u-x*)) 

= g(x*)  + XVg(x*) • (u-x*)  + Xa(^) 
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Here  a:R  -» R is  a function  which  satisfies  lim^  a(X)  = 0.  Since 
g(x*)  < g(x*  + A(u-x*))  for  all  X 5 ^0,1],  Vg(x*) • (u-x*)  + a(X)  > 0 
for  all  X c [0,1].  Since  lim^  ct(X)  =0  we  have  Vg(x*)  • (u-x*)  >0. 
Since  this  is  true  for  all  u € K,  (x*,Vg(x*) ) is  a stationarj’  point 
of  the  pair  (Vg,K). 

Given  a convex  cone  K in  R*^,  its  polar  is  the  set  of  all 

vectors  y in  R*^  such  that  x-y  > 0 for  all  x £ K.  If  f;K  -»  r"^ 
then  the  complementarity  problem  (f,K)  is  to  find  (x,y)  satisfying 

X £ K,  y £ f (x)  n and  x*y  = 0 

Much  effort  has  been  expended  in  studying  the  coII^)lementarity  problem. 
Moreover,  a wide  range  of  applied  problems  have  been  reduced  to  it. 
However,  in  many  instances  the  stationary  point  problem  is  the  more 
convenient  reduction.  The  next  result  from  Saigal  [12]  relates  the 
two  problems. 

n*  / 

Proposition  1.  Suppose  K is  a convex  cone  and  f;K  ->  R . Then  (x,y) 
solves  the  complementarity  problem  if  and  only  if  it  solves  the  stationary 
point  problem. 

Proof.  Suppose  (x,y)  solves  the  stationary  point  problem.  Then 
y £ f(x)  and  x-y  < u-y  for  all  u £ K.  Since  0 £ K,  X'y  < 0.  Also 
X £ K implies  px  € K for  p > 0.  So  for  p > 0,  (p-l)x'y  > 0. 

Since  p-1  changes  sign,  we  must  have  x-y  = 0,  and  thus  u-y  > 0 for 
all  u £ K which  says  that  y £ K^,  so  (x,y)  solves  the  complemen- 
tarity problem. 

Next  let  (x,y)  solve  the  complementarity  problem.  Hence  for 
all  u £ K,  u-y  > 0 = x*y  or  (u-x)-y  > 0.  Now  I need  only  to  note 
that  X £ K and  y € f (x) . □ 
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The  following  proposition  concerns  another  important  special 
case  of  the  stationary  point  problem.  It  also  illustrates  the  relation- 
ship between  the  stationary  point  problem  and  stationary  points  of  mathe- 
matical programming  problems. 

Proposition  2.  Suppose  A £ b £ R*^,  and  K = {x£r'^|Ax  < b] . 

Let  f;R°  R^  . Then  (x,y)  solves  the  stationary  point  problem 
(f,K)  if  and  only  if  x £ K,  y £ f(x),  and  there  is  a v > 0 for  which 
TT-  (Ax-b',  = 0 and  A tt  + y = 0. 

Proof.  Let  (x,y)  solve  the  stationary  point  problem.  Then  x solves 
this  linear  program. 

min  yu 
Au  < b 

Let  ^ be  the  solution  to  the  dual 

max  -b-TT 
T 

A n + y = 0 
7T  > 0 

T 

Now  •n.(Ax-b)  = 0 since  A v + y = 0 and  y-x  = -b^v. 

T 

Conversely  if  Ax  < b,  r > 0,  n'(Ax-b)  = 0,  and  A tt  + y = 0 
then  -b*-"  = -jr-Ax  = y-x  so  x solves  the  primal  linear  program.  If 
in  addition  y £ f(x)  then  f solves  the  stationary  point  problem.  □ 

It 


2.2.  Existence  Theorems  j 

This  section  has  two  main  theorems  which  deal  with  the  stationary 
point  problem.  Theorem  5 is  a very  general  but  non-constructive  result. 

The  proof  is  adapted  from  Saigal  [!’],  but  the  theorem  is  more  general 
in  that  it  deals  with  the  stationary  point  problem  rather  than  the 
complementarity  problem.  If  A and  B are  sets  then  we  write  A'V  B 

i 

for  the  intersection  of  A and  the  complement  of  B.  t 

] 

Theorem  j.  Let  K C be  closed  and  convex,  C C K be  compact,  and 

-X-  ' 

let  f;K  -»  r”  be  upper  semicontinuous.  Suppose  that  for  all 

xc  K\C  , y £ f(x)  there  is  a u 1 C for  which  (x-u)*y  > 0 . | 

Then  the  stationary  point  problem  (f,K)  has  a solution. 

Proof.  Define  D to  be  any  compact  convex  set  C C D c K such  that 
for  all  X € C,  y £ K\  D there  is  a A > 0 such  that  (l-A)x  + Xy  £ D. 

Let  E be  any  compact  convex  set  containing  f(D). 

For  each  y £ r"^  define  iriy)  to  be  the  set  of  ali  x £ D 

which  minimize  X'y  on  D.  Then  -n  is  upper  semicontinuous  and  for 

each  y £ R*^,  m(y)  is  nonempty,  convex,  and  compact. 

Define  PtD  x E -x  D*  x E*  by  (x,y)  -x  {jr{y) , f(x)).  By  the 
Eilenberg-Montgomery  fixed  point  theorem  [ 8 ],  D has  a fixed  point. 

Call  it  (x,y). 

I now  need  only  show  that  x'y  £ x-y  for  all  x Z K.  Suppose 
to  the  contrary  that  there  is  some  x £ K with  x-y  < x-y.  Then 
X K\  D.  Also  X ^ C since  if  x C then  there  is  a X > 0 such 

that  (l-X)x  + Xx  ..  D,  but  this  implies  [ (l-X)x  + Ax]‘y  < x-^  which 

contradicts  the  definition  of  x.  So  x / C. 

lb 


Then  there  is  a u c C with  (x-u)-y  >0  or  x*y  > u-y.  There 
is  a X > 0 with  (l-A)u  + Xx  € D,  so 

* (l-X)u  + AxJ-y  = (l-Aiu’y  + Ax-y  < x-y 

which  contradicts  the  definition  of  x.  So  for  each  x £ K,  x*y  < x*y 
and  (x,y)  is  a stationary  point  of  (f,K).  □ 

The  basis  for  the  constructive  results  of  this  section  is  the 

algorithm  of  Eaves  and  Saigal  [7]  for  systems  of  nonlinear  equations. 

The  next  lemma  gives  a general  condition  on  a point-to-set  map 
n n*  n 

f;R  R under  which  an  x c R may  be  approximated  which  satisfies 
0 £ f(x).  The  algorithm  solves  the  system  of  equations  using  piecewise 
linear  approximations.  The  solution  is  approached  by  solving  finer  and 
finer  approximations. 

We  say  U separates  v ^ U from  infinity  if  every  connected 
set  containing  v and  not  meeting  U is  bounded. 

Let  C C"  R^.  We  define  conv(C)  to  be  the  intersection  of  all 
convex  sets  containing  C.  We  write  dia(C)  for 

sup  ||x-y||  . 
x,y€C 

In  proving  the  following  theorem  we  list  some  of  the  most  important 
properties  of  the  Eaves-Saigal  algorithm.  A complete  description  of 
the  algorithm  may  be  found  in  Eaves  and  Saigal  f7]. 
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Lewmn  Let  ftR*^  ->  r”*  be  upper  semicontinuous.  Let  v £ R*^. 

Let  U - r"^  be  a closed  set  which  separates  v from  infinity.  Let 
D be  the  union  of  all  connected  sets  containing  v but  not  meeting  U. 
Define  G:R^  ->  by  g(x)  = conv(f(x)  U (x-v)).  Suppose  that 

0 £ G(U).  Then  there  is  a & > 0 such  that  if  the  initial  grid  size 
is  less  than  5 then  the  Eaves-Saigal  algorithm  will  not  evaluate  f 
outside  D.  Therefore  it  will  converge  to  a connected  set  f 0 
with  0 f(x)  for  all  x £ C . 

Proof.  Given  a starting  point  v and  an  initial  grid  size  5 the 
Eaves-Saigal  algorithm  generates  a sequence  of  sets  i = 0,1,2,... 

with  the  following  properties. 

(2.1)  V £ conv(CQ) 

(2.2)  dia(C^)  < 6 

(2.3)  0 £ conv(f(C^)  U (C^-v)) 

(2.U)  = C.  U {s^)\{t.] 

(2.5)  if  C ^ X on  a subsequence  then  0 £ conv(f(C.  ))  for  all 
large  j . 

(2.6)  If  all  lie  in  a bounded  set,  then  dia(C^)  -» 0. 

Let  C = Tim  C,.  If  x G C then  0 £ f (x) . If  all  the  C. 
00  1 00  1 

lie  in  a bounded  set  then  is  a nonempty,  closed,  connected  set. 


• “ 1 .-A  j.  * •Ctari.C'*:  • 
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D is  a connected  set  containing  v,  and  not  meeting  U,  so  D 
is  bounded  since  U separates  v from  infinity. 

We  now  use  an  indirect  argument  to  verify  the  theorem.  If  the 
conclusion  does  not  hold  then  there  is  a sequence  6 -»  0 such  that 

for  each  k the  algorithm  produces  a sequence  i = 0,1,2,...  which 

is  not  contained  in  D.  For  each  k let  x be  the  first  point  in 

the  sequence  C‘^  which  is  not  in  D.  By  (2.4)  - {x  ] c D so 

d(x^,D)  < 5^. 

*"k  k ***  k k k 

Let  X be  any  point  in  D with  d(x  ,x  ) < & . Since  x ^ D 

there  is  a point  u^  € conv(x*'’,x^)  H U.  Therefore  d(x^,U)  < 5^. 

Let  X*  be  any  cluster  point  of  the  x . Then  d(x*,U)  = 0, 
so  X**'  1 U since  U is  closed.  But  since  for  all  i and  all  k, 

0 £ conv(G(C^)),  it  follows  that  0 £ G(x*).  This  contradiction  proves 
the  theorem.  □ 

In  practice  neither  6 nor  the  set  U will  be  known  to  the  user 
of  the  algorithm.  However,  if  it  is  knovm  that  the  function  f and 
starting  point  v satisfy  the  hypotheses  of  the  lemma  for  some  U 
then  the  following  process  may  be  employed  to  find  a solution. 

1. ^  Choose  and  M^,  k = 0,1,2,...  such  that  0 and 

M Set  k = 0. 

2. )  Apply  the  algorithm  to  f using  starting  point  v and  initial 

k ri 

grid  size  6 . If  no  points  x £ R are  generated  with 

||x-vl|  > then  the  algorithm  must  converge  to  a set  of  solutions. 

OtheiT/ise  go  to  3. 
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5.'  Set  k = k+1.  Go  back  to  2. 

Eventually,  a k will  be  found  for  which  6 is  small  enough 
and  l!x-vll  < for  all  x C D.  Lemma  U tells  us  that  the  algorithm 
will  converge  to  a set  of  solutions  for  this  k. 

In  practice  grid  size  does  not  appear  to  be  critical.  If  the 
algorithm  seems  to  be  straying  too  far  from  the  starting  point,  it  would 
probably  be  better  to  look  for  errors  in  the  formulation  than  to  try 
again  with  a smaller  grid  size. 


Lemma  s.  If  p,  q,  and  r are  in  and  ||p-q||  < ||p-r||  then 

';p-r)(q-r'  >0. 

Proof.  !lp-q||^  < llp-r||^ 

ll(p-r)  - (q-r)ll^  < ||p-r||^ 

l|p-r||^  - ?(p-r)(q-r)  + ||q-r||^  < ||p-r|i^ 

-2(p-r)  (q-'r'  + ||q-r||^  <r  0 

(p-r)(q-r)  > 0 . □ 

The  next  theorem  is  a constructive  version  of  Theorem  3.  In 
order  to  be  able  to  compute  a solution  we  need  only  slightly  stronger 
conditions  than  in  Theorem  3. 

A function  grR^  -»  R is  convex  by  definition  if  for  all  x,  y £ R*^, 

- ^0,1],  g(Ax  + (l-X)y)  < Xg(x)  + (1-X)  g(y).  If  g is  convex  then 
we  define  the  subgradient  map  of  g by 
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'^g(x)  = fy  R'^ly  (x-u)  > gfx)  - g(ui  for  all  u € r'^]  . 


If  g is  convex  then  dg(x)  is  nonempty,  compact  and  convex  for  each  x 
and  'g  is  upper  semicontinuous. 

Theorem  b.  Let  -»  R i = 1,  . . . , ra  be  convex.  Let 

K = {x  e R'^lg^(x)  <0,  i = l,...,m].  Let  C C K be  compact.  Let 

n 0 

f;R  -^R  be  upper  semicontinuous.  Suppose  that  there  is  a u £ C 

such  that  gj^(u*^)  <0,  i = 1,  Suppose  that  for  all 

X £ K\C,  y £ f(x)  there  is  a u £ C such  that  (x-u'  -y  > 0.  Then 

the  stationary  point  problem  (f,K'  has  a solution  which  can  be 

computed  by  the  Eaves-Saigal  algorithm  using  any  starting  point. 

Proof.  For  each  x £ R*^  define  F(x)  to  be  the  set  of  all 
Z'!'  „ A.y^  which  satisfy  these  conditions. 

m 

Z A.  = 1 A > 0 

i=0  ^ 

y°  £ f(x)  ^0  ^ ° i ^ 0,  gj^(x)  > 0 

(2.7)  For  i = 1, . . ,m 

y^  £ ag^(x)  Aj^  = 0 if  g^(x)  < 0. 

First  we  need  to  show  that  F is  upper  semicontinuous.  Suppose 
x"^  -»  x,  and  w"^  € F(x^)  for  each  j.  For  each  j there  exist 


that 


and  i = which  satisfy  (2.7)  and  such  that 

Since  the  converge,  they  lie  in  some  compact 

set.  The  y^’^  are  in  the  (compact)  images  of  this  compact  set 
under  f,  )g^,  ...  , Thus  there  is  a subsequence  on  which  the 

and  converge,  say  to  y^  and  A.  In  the  following  we 

will  deal  only  with  this  subsequence. 

Since  the  are  continuous,  for  all  .j  sufficiently  large 


and  for  all  i = 1, . . . ,m,  g 

. (x^)  >0 

X 

if  g^(x)  > 

0,  and  g^( 

V 

X 

if  g^(x)  < 0.  If  for  any 

H* 

II 

,m,  g^(x)  > 0 

then  Aq 

= 0 : 

large  j,  so  Aq  = 0.  If 

II 

m and  g^(x) 

< 0 then 

A^  = 
1 

for  large  j,  so  A.  = 0. 

J 

Since  f 

and  the  Sg^^ 

are  upper 

semi- 

continuous,  y°  € f(x)  and  y^  € <5gj_(x).  Therefore  A,  x,  and  y^, 
i = 0,...,m  satisfy  (2.7).  Let  w = Then  w£  F(x)  and 

on  the  subsequence  w"^  w so  F is  upper  semi  continuous. 

It  is  clear  that  F(x)  is  nonempty,  convex  and  compact  for 

each  X. 

Let  V Let  M = sup^^  l|u-v||.  To  show  that  F satisfies 

the  hypotheses  of  Theorem  4 it  only  remains  to  show  that  for  all  x € 
with  l|x-v||  > M that  0 € conv(F(x)  {x-v}). 

So  suppose  X £ r”^  and  |!x-vil  > M.  Let  w £ F(x'  . There  exist 
A,  y^...-,y'”  which  satisfy  (2.7)  such  that  w = E"?  „ A.y^.  There  is 
a u^  £ C such  that  g^(u^)  < 0,  i = 1, ...,m.  We  distinguish  two 
cases  according  to  whether  = 0. 

If  Aq>0  then  x£K\C  so  there  is  a uIC  such  that  (x-u)-y*^  > 0. 
Choose  a sufficiently  small  that  (x-(cai^  + (l-a)u))-y^  > 0.  Let 
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u'  = ou®  + (l-a)u.  Then  g^(u' ) < 0,  i = 1,  . . . ,ni.  For  each  i = 1, . . . ,m 
if  > 0 then  g^(x}  > 0,  so  bv  the  gradient  inequality  (x-u')'y^ 

> g^(x)  - g^(u''i  > 0.  Now  (x-u'  )-w  = A^(x-u')-y^.  This  sum  is 

positive  since  no  term  is  negative  and  some  term  is  positive.  Since 
j|x-v|l  > |)v-u'||  we  have  (x-u')-(x-v'  >0  by  Lemma  5-  Therefore 
0 ^ conv((w]  U (x-v)). 

But  if  Aq  = 0 then  for  each  i = 1,  ...,m  if  > 0 then 
g^(x^  >0  so  by  the  gradient  inequality  (x-u*^)-y^  > g^(x)  - g^(u^)  > 0. 
Now  (x-u^l'W  = A^(x-u*^)y^.  This  sum  is  positive  since  some  term 

is  positive  and  no  term  is  negative.  Since  ||x-v||  > liv-u^||  we  have 
(x-u*^)  • (x-v)  >0  by  Lemma  5.  Therefore  0 i/  conv((x)  U {x-v}). 

Since  0 ^ conv({w]  U {x-v]  ) for  all  w € F(x'j, 

0 s'  conv(F(x)  U (x-v}). 

By  Theorem  U we  may  use  the  Eaves-Saigal  algorithm  to  approximate 
an  X € with  0 € F(x).  For  such  an  x there  exist  A and  y^, 

1 = 0,...,m  which  satisfy  (2.7)  and  ^_q  A^y^  = 0.  For  each 

i = 1,  ...,m  if  > 0 then  g^(x'  > 0 so  by  the  gradient  inequality 
(x-u'^)y^  > gi^(x)  - g^(u*^)  >0.  If  Aq  < 1 then 

A^(x-u°)y°  = - Z A. (x-u°)y^  < Z A. (g. (u°)  - g. (x) ) < 0 . 

i=l  ""  i=l  ^ 

Thus  ^ u € K.  Then 

Ar.(x-u)y^  = - Z A.  (x-u)y  < Z A.  (g,  (u)  - g.  (x) ) < 0 . 

i=l  i=l 

So  (x-u)y^  <0  and  (x,y^)  is  a stationary  point  of  (f,K^.  □ 


.3.  Uniqueness  Theorems 


In  this  section  we  explore  the  question  of  unioueness  of 
stationary  points. 

n 

Let  f;K  ->  R where  K c R is  closed  and  convex.  The 
mapping  f is  called  monotone  on  K if  for  every  pair  of  vectors  x" 
and  x^  in  K and  every  y^  f f(x^)  and  y^  £ f(x^i,  (x^-x'^)  • (y^-y^)  > 0. 

The  mapping  f is  strictly  monotone  if  for  every  pair  of  distinct  vectors 

IP  11  2^2 

X and  x“  in  K,  and  every  y G f(x  ) and  y £ f(x  ), 

IP  12 

(x  -X  ''•(y  -y  ) >0.  The  following  is  a useful  uniqueness  result  for 
the  stationary  point  problem. 


THEOREM  7.  Let  f:K  -4  R be  strictly  monotone.  If  there  is  a solution 
to  the  stationary  point  problem  (f, K)  then  it  is  unique. 

"^1  p 2 11 

Proof.  Suppose  ) and  (x‘',y  ) solve  (f,Ki.  Then  y £ f(x  ), 

y^  £ f(x^'  and  for  all  u £ K,  (u-x^)-y^  > 0 and  (u-x^i -y^  > 0. 

2 11  12  2 

In  particular  (x  -x  )-y  > 0 and  (x  -x  ).y  > 0.  Summing  these 

1 P 1 2 

we  have  (x  -x  ) (y  -y  ) < 0.  But  this  is  impossible  since  f is 
strictly  monotone.  □ 

The  following  theorem  is  based  on  the  same  idea,  but  is  slightly 
more  general.  By  using  linearity  and  relaxing  the  strict  monotonicity 
we  obtain  not  a unique  solution,  but  a convex  set  of  solutions.  However, 
in  many  cases  we  may  be  able  to  determine  the  entire  set  of  solutions. 
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THEOREM  8.  Let  K c R*^  and  B £ r”  ^ Let  L = {Bx]x  £ K] . Let 
fjL  r"'  be  strictly  monotone.  Let  A € be  positive  semi- 

definite.  Define  g:K  -»  by  g(x)  = B'^f(Bx)  + Ax.  Then  the  set 

of  solutions  (x,y)  to  the  stationary  point  problem  (g,K)  is  convex 
and  Bx  has  the  same  value  for  all  of  these  solutions. 

11  2 2 

Proof.  Suppose  (x  ,y  ) and  (x  ,y  ) solve  (g,K).  Then  there  exist 

w*  and  such  that  for  i = 1,2,  y^  = b'^w^  + Ax^,  € f(Bx^)  and 

for  all  u £ K,  (u-x^)  • (b'^w^  + Ax^)  > 0.  In  particular  (x^-x^)  • (B^w^  + Ax^)  >0 

and  (x^-x^)  • (b'^w^  + Ax^)  >0.  Summing  these  (x^-x^)  • (b'^(w^-w^)+A(x^-x^)  ) <0 

or  (x^-x^) -ACx^-x^)  + B(x^-x^) • (w^-w^)  < 0.  But  since  A is  positive 
12  12 

serai-definite  B(x  -x  )‘(w  -w  ) <0.  Since  f is  strictly  monotone 

1 2 

this  implies  that  Bx  = Bx  , 

Let  A £ [0,1],  X = Ax^  + (l-A)x^,  y = Ay^  + (l-A)y^,  and 

12  1 2 T 1 

w = Aw  + (l-A)w  . Then  Bx  = Bx  = Bx  and  y = B w + Ax,  Since  w 

p 

and  w are  in  f(Bx)  and  f(Bx)  is  convex,  w € f(Bx).  So  y £ g(x). 

Let  u £ K.  Since  A is  positive  semi-definite 

(u-x)Ax  > a(u-x^)Ax^  + (l-A) (u-x^)Ax^ 

(u-x)b'^w  = A(u-x)B^w^  + (l-A)  (u-x)b'^w^ 

= A(u-x^)B^w^  + (1-A)  (u-x^)b'^w^ 

(u-x)(b'^w  + Ax)  > A(u-x^)  (b'^w^  + Ax^)  + (l-A)  (u-x^Hb'^w^  + Ax^)  > 0 
solves  (g,K)  and  the  theorem  follows.  □ 


pu 


Therefore  (x,y) 
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PROOFS  OF  THE  ^1AIN  THEOREMS 


This  chapter  is  devoted  to  the  proofs  of  the  theorems  of 
Chapter  L.  Theorem  1.1  will  not  be  given  a separate  proof  since  it 
is  clearly  a special  case  of  Theorem  1.2.  Here  t > e means  no 
component  of  t is  less  than  6 . 


Proof  of  Theorem  1.2. 

Throughout  this  proof  tj^  ^ will  represent  the  constsint  zero. 
Thus  t >•“  where  m = |^  |•('|^|-l).  For  computational  efficiency 
we  will  make  the  additional  requirement  that  y. . . =0  for  all  arcs 
ij.  In  the  following  we  will  show  that  a solution  to  the  equilibrium 
problem  can  be  computed  which  satisfies  this  restriction.  For  this 
reason  we  will  regard  y as  a vector  in  r'^  where  n = q*  |-1' 
and  q = |fl|  is  the  total  number  of  arcs.  Define  B t R ^ ^ 


by 


ij  € a . 


Then  we  may  write 


(3  x) . . = X.  . 

ij 


i ^ K ij  e a,  k e 


Define  a matrix  A 6 by 


i ^ K,  i,k  e 


2^ 


Then  we  may  write 


^i,k  “ ^j,k 


i ^ k,  ij  f-:  fl,  k £ ^. 


In  these  formulas  y. . . and  t are  to  be  read  as  zero 

1J,1  K,K 


where  they  appear. 


For  each  i ^ k,  £,  k t choose  a path  from  £ to  k.  Let 

k n k 

d ’ be  a vector  in  R such  that  d.’.  . = 1 if  ij  is  on  the  path 

ij  j k 

£ k i k 

and  all  other  components  of  d ’ are  zero.  Thus  Ad  ’ is  a vector 

in  r”  such  that  (Ad”^’^)  =1  and  all  other  components  are  zero. 

K , K 

Choose  b £ R™  so  that  b > g(t)  for  all  t > c.  Choose 

h £ r”  so  that  for  all  ij  £ 3,  k £ 0^,  i ^ k 


Let  K = {(y,t)  £ R^'^’^jt  > e,  0 < y < h] . Define 


_ „m+n  „m+-n*  , 

G:  R R by 


G(y,t)  = 


'B'^f(By)  -A'^t 


Ay  - g(t) 


We  proceed  to  show  that  the  stationary  point  problem  (G,K) 
may  be  solved  as  described  in  Chapter  2,  and  that  this  solution  will 
also  solve  the  equilibrium  problem. 

Let  y^  = it  Then  for  each  ij  £ 3, 

i ^ k - ^ we  have  the  following. 
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'iO,k.  ^ ' £,k  -..j,-  ^ 

y®  < h 


Ay°  = Z (b.  . + l)Ad^’^  = b + e 
i,k 


Here  e € r”  is  a vector  of  ones.  Let 


M = -inf{(y-y°)-B’’T|0  < y < h,  x £ f(By)) 


Then  M > 0 since  0 < y < h.  Let 


C = {(y,t)  G R'^^'Ky.t)  € K,  e-t  < M} 


If  (y,t)  e k\^C  then  for  each  r ^ g{b) , x £ f(By) 


y - y 


T T 
■^T  - A t 


Ay  - r 


(y-y°)-B^T  + t-(Ay°-r) 


> -M  + e-t  > 0 . 


So  by  Theorem  2.6  the  stationary  point  problem  (G,K)  can  be 


solved. 


Suppose  ((y,t) , (B'^T,r ' ) solves  the  stationary  point  problem 


(G,K).  By  Proposition  2.2  we  may  write  this  in  the  following  way. 

t>e,  0<y^h,  rt  g(t),  T ' f(By) 


’ i,k  - ^ ^ij,k  ■ j ^ji.k 
^^i.k''^'^‘^^i,k  ■ ^ yij,k  ^ ^Ji.k'  = ° 


i k,  i.  k e ^ 
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i j € a,  k t i ^ k 


[ . . + t - t 
ij  i,k 


>0  if  y.  . ^ < h.  . ^ 
ij,k  ij,k 


<0  if  y.  . , > 0 


If  for  some  l,k,  - Ej  y_,i  then  t.  = e 

and  for  all  arcs  in  in  along  some  path  from  i to  k,  y.  . , >0. 

1 d ^1^2’^ 

Therefore  t.  , > t.  , +€  since  r.  . > • . But  since  t,  =0, 

i^,k  i2,k  1^’-  k,k 

■t.  v > ty  transitivity  along  the  path.  So  r.  , = S.  y.  . , - Z.y..  , 

1’*^  i,k  j ij,k  j ji,k 


after  all. 


If  for  some  i j c a,  k £ ^ , i k,  y . = h 

• 1 J ^ JV  1 1]  > J 


then 


^ii  k ^ k - k'  traffic  on  arc  ij  with 

destination  k is  greater  than  the  traffic  from  all  origins  with 

destination  k.  This  is  only  possible  if  there  is  a loop  which  includes 

ij  such  that  for  all  i i in  this  loop  y.  > 0-  This  implies 

d ^1^2’^ 

that  t.  , > t , since  t . >0.  But  this  is  impossible  since 

Igyk  i^i^ 

it  implies  by  transitivity  along  the  loop  that  t.  , is  less  than 

1,  K 

itself.  Therefore  y < h. 

Since  t > e it  follows  that  t > 0.  Since  t > 0,  t > 0 

and  t,  , = 0 for  all  k,  t.  + t - t > 0. 
k,k  ij  J,i  1,1 

Taking  the  last  three  paragraphs  into  account  we  may  write  the 
solution  conditions  for  the  stationary  point  problem  in  this  ways 


t > 0, 


y > 0,  re  g(t), 


r e f(By) 


'^i,k  = ^ ^ij.k  " J yji,k 


2P 


i ^ k,  i,  k £ ^ 


ij  e a,  k L 


r.  . + 
ij 


r.  . 

ij 


+ V 


‘l,k 


But  these  are  precisely  the  desired  equilibrium  conditions  if  we  let 
X = By. 


Before  proving  Theorem  1.3  we  will  need  the  following  lemma 
on  the  minimum  of  a continuous  function.  If  fiR^^  -»  R and  D C R^ 
then  we  define 


arg  min  f(x)  = {x  £ D|f(x)  < f(y)  for  all  y £ D] 
x£  D 


ni  * 

Lemma.  Let  K:R  -> D be  continuous,  where  D is  a compact  subset 

of  R*^.  Let  f ; R*^  X R^  ->  R be  a continuous  function.  For  each 

z £ R*”,  w £ R^  let  j(z,w)  = arg  min  f(x,w).  Then  J is  upper  semi 

x£K( z) 

continuous. 

Proof.  Suppose  z^  z,  w^  w,  and  x^  £ J(z^,w^)  for  all  i.  Since 
each  x^  £ D there  is  a subsequence  on  which  the  x^  converge,  say  to 
X.  Since  K is  upper  semicontinuous,  x £ K(z).  Hereafter  we  will  work 
only  on  this  subsequence.  If  x £ K(z)  then  since  K is  lower  semi- 
continuous there  is  a sequence  x^  ^ x such  that  x^  £ K(z^)  for  all  i. 
Therefore 

f(x,w)  = lim  f(x^,w^)  < lim  f(x^,w^)  = f(x,w). 
i ~ i 

So  X C j(z,w)  which  implies  that  J is  upper  semicontinuous.  □ 


Proof  of  Theorem  1.3.  For  each  j = and  each  t > e let 


= {(q^r'^) \r^  >0,  < q^  < q^}  . 

Then  for  each  (q'^,r^)  £ K^(t)  we  have  0 < < q^  and 

e-r"^  < - t-y'^  < i < i q'^'  , 

— e — € — € 

Therefore  U K‘^(t)  is  bounded.  Define 
t >e 

h'^(t)  = arg  min  U(q'^, 

(q^y'^)£K'^(t' 

Then  n'^  is  bounded  on  [ 1 1 1 > e } . 

For  each  t > e there  exists  (q'^,y'^)  £ K'^(t)  such  that  y*^  > 0, 
q"^  < q'^  and  y^-t  < q^.  So  by  [6,  Corollary  II. 3. 2]  is  continuous. 

By  the  lemma  h^  is  upper  semicontinuous . But  g*^  can  be  obtained  by 
deleting  the  first  component  in  the  range  of  h'^.  Therefore  g*^  is 
nonnegative,  bounded,  convex  valued,  and  upper  semicontinuous  on 
(t|t  > f],  and  so  is  g since  g = Z.g^.  □ 

Proof,  of  Theorem  1.4.  In  the  proof  of  Theorem  1.2  we  introduced  a 
stationary  point  problem  (G,K)  and  showed  that  solutions  to  (G, K) 
are  also  solutions  to  the  network  equilibrium  problem.  Here  we  show 
the  opposite,  namely  that  solutions  to  the  equilibrium  problem  solve 
the  stationary  point  problem.  The  theorem  will  follow  when  we  have 
shown  that  (G,K)  satisfies  the  hypotheses  of  Theorem  2.8. 
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Since  g is  nonnegative  and  -g  is  strictly  monotone,  g 
is  positive.  Suppose  y,  t,  x,  and  r solve  the  equilibrium  problem. 
We  may  state  this  as  follows. 

t > 0,  y > 0,  r £ g(t),  T £ f(By) 


^i,k  “ ^ ^'ij,k  ■ ^ yji,k 

<J  t) 


i k,  i,  k £ 7^ 


T..+t.,-t.  ,>0 
10  o,k  i,k  - 


(T..  + t..-t.  ,,)-y.  = 0 

10  0,k  i,k  ■'io,k 


ij  £ a,  k 


For  any  pair  of  nodes  i,  k since  u ^ 0 then  for  all 

1,  K 

arcs  i, io  along  some  path  from  i to  k,  y . ^ > 0.  Therefore 

\,k>%,k'"^  V2^"'  \,k=°’^,k>" 

by  transitivity  along  the  path.  So  t > e. 

If  for  some  i j £ fl,  k £ ^ , i ^ k,  y, . , > h ^ then 
^ij  k ^ k - k"  traffic  on  arc  ij  with 

destination  k is  greater  than  the  traffic  from  all  origins  with 
destination  k.  This  is  only  possible  if  there  is  a loop  which  includes 
ij  such  that  for  all  1,1^  in  this  loop  y,  . . > 0.  This  implies 

1 1 J ^ 

that  t 1,  > t ^ since  x.  . >0.  But  this  is  impossible  since 

it  implies  by  transitivity  along  the  loop  that  t , is  less  than 

1,  K 

itself.  Therefore  y < h. 

Since  t > € and  y < h,  we  have  in  fact  a solution  to 

(G,K). 
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To  show  that  the  hypotheses  of  Theorem  2.8  are  satisfied  we 


make  the  following  definitions. 


f(x,t)  = (f(x),  -g(t) ) . 

-T-  — — — 

Clearly  G(y,t)  = B f(B(y,t))  + A(y,t).  f is  strictly  monotone  because 
f and  -g  are  strictly  monotone.  It  is  easy  to  see  that  A is 
positive  semidef inite.  Therefore  the  set  of  solutions  to  the  stationary 
point  problem  (and  the  equilibrium  problem)  is  convex,  and  B(y,t) 
has  the  same  value  for  all  such  solutions.  In  particular  t is  unique 
Sind  By  is  unique  but  the  flows  y may  not  be  unique.  □ 


CHAPTER  i; 


COMPUTATION 

The  purpose  of  this  chapter  is  to  demonstrate  the  computational 
viability  of  the  procedures  described  in  Chapters  2 and  3, 

^4. 1.  Formulation  and  Approach 

The  method  used  to  solve  the  network  problems  is  outlined 
in  the  proofs  of  Theorems  I.  , 2.h  and  2.6.  The  proof  of  Theorem  1.2 
(Chapter  3)  shows  how  the  network  problem  can  be  formulated  as  a 
stationary  point  problem  which  satisfies  the  conditions  of  Theorem  2.6. 

The  prcof  of  Theorem  2.6  shows  how  such  a stationary  point  problem 
can  be  formulated  as  a problem  of  finding  a zero  of  a point-to-set  map. 

H.  Z.  Aashtiani  [1]  has  had  some  success  with  an  approach  to 
solving  traffic  network  problems  which  uses  Lemke's  algorithm  as  a sub- 
routine. However,  no  convergence  proof  is  offered  for  this  process. 

Considerable  latitude  is  available  in  solving  a stationary  problem 
(f,K'  using  the  Eaves-Saigal  algorithm.  The  remainder  of  this  section 
is  devoted  to  a trick  vrhich  uses  this  latitude  to  reduce  the  time  required 
to  solve  certain  stationary  point  problems. 

In  the  proof  of  Theorem  2.6  we  defined  a map  F by  setting 

HI 

Ffx)  equal  to  the  set  of  all  such  that 

m 

E A,  = 1 A > 0 

i=0 

y°  e f(x)  ® > 0 
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For  i = 1, . . . ,m 


t ^g^(x)  = 0 if  g^(x)  < 0. 

In  solving  a point-to-set  map  F using  the  Eaves-Saigal  algorithm 
a unique  member  of  F(x)  must  be  chosen  for  each  point  x t It 

is  often  the  case  that  f and  c'g^  are  continuous,  so  that  the 
are  unique.  However  we  still  must  choose  the  weights 

In  addition  we  may  choose  the  triangulation  or  grid  on  which 
the  Eaves-Saigal  algorithm  operates.  In  the  case  where  the  functions 
g^  are  affine  it  is  possible  to  choose  a triangulation  in  which  every 
simplex  which  meets  the  interior  of  K is  contained  in  K.  Using  such 
a triangulation  with  a choice  of  A where  Aq  = 1 if  x £ K,  we  obtain 
greatly  superior  convergence.  This  is  due  to  the  fact  that  a solution 
to  the  piecewise-linear  approximation  at  large  grid  sizes  is  a good 
approximation  to  the  actual  solution.  This  is  not  true  for  most  choices 
of  weights  A and  of  the  triangulation.  Finding  a suitable  triangulation 
is  particularly  easy  in  the  case  where  the  constraints  on  K are  only 
upper  and  lower  bounds,  as  they  are  for  the  network  problems  of  Chapter  1. 

In  the  next  section  we  will  demonstrate  empirically  the  usefulness 
of  this  technique  in  the  case  of  the  traffic  flow  problems. 

h.2.  Examples. 

In  this  section  two  examples  are  given  of  the  network  formulation 
of  Chapter  1.  In  order  to  solve  these  problems  we  prepared  a computer 
code  which  provides  an  interface  between  the  network  problem  and  the 
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Eaves-Salgal  algorithm.  We  used  a program  by  Romesh  Saigal  to  execute 
the  Eaves-Saigal  algorithm. 

All  of  the  computer  runs  were  made  on  an  IBM  3?0/l68  using 

FORfRAN  H with  the  highest  level  of  optimization. 

If  D C R^  then  the  diameter  of  D is  sup  T^llx-y||-  The 

x,y_D 

initial  grid  sizes  reported  below  refer  to  the  diameter  of  the  simplices 
on  ’-diich  the  function  is  first  evaluated.  The  final  grid  sizes  refer 
to  the  diameter  of  the  final  simplex. 

The  first  example  is  based  on  the  following  network. 


Here  = (1,2,5)  and  a = (12,21,25,51) . The  delay  fianctions  are 


fl2(y)  = 10  + e 


+ 1.25  log(y^^  + 1.0) 


where 


fp^Cy)  = ^ + e 
+ e 


The  travel  demand  functions  are  as  follows. 


80 


j + 1 


*1,3'’^' 


120 


Uo 


gg^iCt) 


Vi""- 


100 


^2, 1 — ^2,3 


if  t^  - < t_  , 

2,1  - 2,3 


80 


«2,3<*> 


‘2,3^^ 


20 


‘2,3*^ 


if  ^2,1  -^2,3 


V1-V3 


60 


83  1^")  = 

^,i  ^3,1+ i 


100 


Where  more  than  one  function  value  is  given  at  a particular 
point  (e.g.  t_  , = i'Q  value  of  g is  the  convex  hull  of  the 

two  values.  In  this  case  some  of  the  travelers  from  node  2 will  go  to 
either  1 or  3 depending  on  which  is  closest.  If  the  travel  times  are 
equal  then  those  travelers  who  want  to  go  to  either  1 or  3 will  be 
divided  between  the  two  destinations. 


•»  AMW 
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Arcs  12  and  21  represent  a two-way  street.  Note  that  although 


the  function  is  symmetric  with  respect  to  the  two  arcs,  the  Jacobian 
matrix  will  not  be  symmetric  unless  the  flows  in  the  two  directions  are 
equal.  Therefore,  the  method  proposed  by  Dafermos  [it]  would  not  be 
applicable  here.  Furthermore  the  demand  functions  are  neither  invertible 
nor  continuous.  Therefore  the  full  power  of  Theorem  1.1  is  required  to 
solve  this  problem. 

Here  are  the  equilibrium  values. 


i,k 

1,2 

19.50 

5.94 

1,5 

28.45 

4.08 

2,1 

15.22 

2.81 

2,5 

9.15 

7.^ 

5,1 

4.09 

11.79 

5,2 

25.58 

4.10 

f.  . 

i-J 

12,2 

8.04 

19.50 

12,5 

4.08 

19.50 

21,1 

1.15 

15.22 

21,5 

0.00 

15.22 

25,1 

1. 66 

9.15 

25,5 

11.97 

9.15 

51,1 

15.i+6 

4.09 

51,2 

4.10 

4.09 
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This  problem  was  nm  with  and  without  the  special  alignment 
discussed  in  Section  1.  In  each  run  the  initial  grid  size  was  l8.7 

_7 

and  the  final  grid  size  was  2.8  x 10  . The  sajae  answer  was  found  in 

each  case. 

In  the  first  run  the  triangulation  was  aligned  with  K.  In 
this  case  the  algorithm  required  886  function  evaluations,  979  pivots,  and 
1.70  seconds  CPU  time.  The  starting  vector  had  1.25  in  each  component. 

In  the  second  run  no  alignment  was  used.  This  time  the  algorithm 
required  1598  function  evaluations,  16qo  pivots  and  2.83  seconds  CPU 
time.  The  starting  vector  had  1.00  in  each  component. 

The  second  example  is  based  on  the  following  network. 


Now  9/  = {1,2,3)  and  fl  = {12,21,23,32,31).  The  delay- 
functions  are 


■12 


21 


23 


32 


'51 


10  + e 


(,12-10) 


+ 1.25-log(y2^  + 1) 


(yp,-io) 

10  + e + 1.25'log(y^2 


5 + 0.5-e 


(ypj-io) 


+ 0.625- log(yj2 


(y,p-io) 

5 + 0.5-e  + 0.625- log(y, 


23 


15.0  + 1.5e 


(y3i^-5.o) 


+ 1) 
+ 1) 
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where 


' I ''id.k 


The  travel  demand  functions  are  identical  to  the  earlier  example. 
Here  are  the  equilihrixim  values. 


i,k 

^i.k 

1,2 

13.30 

5.59 

1,3 

23.55 

4.89 

2,1 

13.05 

2.85 

2,3 

10.25 

7.11 

3,1 

15.39 

3. 66 

3,2 

8.04 

11. o6 

ij,k 

f.  . 

12,2 

5.59 

13.30 

12,3 

4.89 

13.30 

21,1 

2.85 

13.05 

21,3 

0.00 

13.05 

23,1 

0.00 

10.25 

23,3 

12.00 

10.25 

32,1 

0.00 

8.04 

32,3 

11.06 

8.04 

31,1 

3.66 

15.39 

31,2 

0.00 

15.39 
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This  problem  was  also  run  with  and  without  the  special  align- 
ment discussed  in  Section  1.  In  each  run  the  initial  grid  size  was 

-7 

20.0  and  the  final  grid  size  was  2.9  x 10  . The  same  answer  was  found 

in  each  case. 

In  the  first  run  the  triangulation  was  aligned  with  K.  In 
this  case  the  algorithm  required  1107  function  evaluations,  120U 
pivots  and  2.59  seconds  CPU  time.  Without  alignment  the  algorithm 
required  1921  fxinction  evaluations,  2016  pivots  and  5.92  seconds  CPU 
time.  The  first  starting  vector  had  1.25  in  each  component,  the  second 

1.00  in  each  component. 

From  these  runs  we  may  conclude  that  the  procedure  developed 
here  is  a viable  approach  to  the  traffic  equilibrium  but  that  the 
alignment  technique  significantly  reduces  the  computation  time. 
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